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Major  Goals:  The  objective  of  our  research  is  to  develop  a  set  of  mathematical  tools  that  can  help  solve  design 
and  analysis 

issues  for  dynamics  of  large-scale  interconnected  system  of 

systems.  These  tools  are  a  combination  of  operator-theoretic,  geometric  dynamical 

systems  methods  and  graph-theoretic  methods.  We  extend  the 

analysis  methods  based  on  Koopman  operator  properties  to  include  stable  and  unstable  foliations  and  develop 
methods  in  Generalized  Laplace  Analysis  to  do  so  for  a  large  class  of  nonlinear  interconnected  dynamical 
systems. 

We  will  couple  these  operator-theoretic  methods  to  graph  theoretic  analysis  that  will  rely  on  spectral  analysis 
of  the  interconnection  graph  associated  with  a  dynamical  system  at  hand.  The  coupling  is  enabled  by  the  fact 
that  both  approaches  are  spectral.  The  new  theory  is  able  to  treat  discrete  and  continuous,  stochastic  and 
deterministic,  and  hybrid  systems  within  a  single  framework.  We  pursue  a  research  strategy  to  develop 
effective  algorithms  for  energy  efficiency  in  such  complex  dynamical  systems 
and  test  them  on  an  existing  testbed  site  at  the  University  of  California, 

Santa  Barbara,  that  is  fully  equipped  with  sensors  and  enables  us  to  pursue  detailed  analyses  of  the  real-world 
effectiveness  of  such  algorithms. 

Accomplishments  Under 

Accomplishments:  1 )  Comparison  of  different  dynamical  systems  or  a  dynamical  system  with  data  Is  one  of  the 
core  Issues  in 

both  dynamical  systems  and  control  theory.  The  fruitful  approaches  are  particularly  hard  to  come  by  for 
systems  and  data  that  show  nonlinear  behavior  and  non-Gaussian  noise  characteristics.  We  developed  a 
theory  that  utilizes  spectral  theory  of  linear  operators  (composition,  or  Koopman)  to  provide  the  methodology, 
that  was  originally  formulated  for  measure-preserving  systems.  We  extended  it  to  capture  dissipative  and 
finite-time  dynamics.  The  approach  combines  a  version  of  ergodic  partition  theory  with  Hardy  space  theory  in 
observable  space  (rather  than  in  time). 

2)  We  applied  the  operator-theoretic  viewpoint  to  a  class  of  non-smooth  dynamical  systems  that  are  exposed 
to  event-triggered  state  resets.  The  considered  benchmark  problem  is  that  of  a  pendulum  which  receives  a 
downward  kick  under  certain  fixed  angles.  The  pendulum  is  modeled  as  a  hybrid  automaton  and  Is  analyzed 
from  both  a  geometric  perspective  and  the  formalism  carried  out  by  Koopman  operator  theory.  A  connection  is 
drawn  between  these  two  interpretations  of  a  dynamical  system  by  means  of  establishing  a  link  between  the 
spectral  properties  of  the  Koopman  operator  and  the  geometric  properties  in  the  state-space. 
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3)  We  developed  the  connections  between  the  principle  eigenfunctions  of  the  Koopman  operator  associated 
with  a  diffeomorphism  of  a  Banach  space  with  an  asymp-  totically  stable,  hyperbolic  fixed  point  and  the 
existence  of  a  topological  conjugacy  between  the  diffeomorphism  and  its  linearization  around  the  fixed  point. 

As  long  as  the  principle  eigenvalues  satisfy  the  non-resonance  conditions  of  normal  form  theory,  the  principle 
eigenfunctions  generate  a  commutative  algebra  of  polynomials  from  which  a  sequence  of  coordinate 
transformations  that  linearize  the  diffeomorphism  can  be  generated. 

We  answered  the  following  questions. 

(1)  Is  there  a  sub-algebra  of  C(X)  of  which  every  element  has  an  expansion  into  eigenfunctions  of  Koopman 
operator  associated  with  the  linearization  at  equilibrium? 

(2)  Similarly,  is  there  a  sub-algebra  of  C(X)  of  which  every  element  has  an  expansion  into  eigenfunctions  of 
Koopman  operator  of  the  nonlinear  system? 

(3)  If  the  above  two  sub-algebras  exist,  what  is  their  relation? 

Answering  (1),  the  so-called  principle  eigenfunction  can  be  used  to  generate  an  algebra  of  observables  (called 
the  principle  algebra)  which  have  expansions  in  eigenfunctions.  This  algebra  can  be  pulled  back  to  an  algebra 
of  observables  for  the  nonlinear  dynamics  as  long  as  a  topological  conjugacy  exists.  This  conjugacy  can 
actually  be  constructed  from  change  of  variables  transformations  that  are  formed  from  tensor  products  of  the 
principle  algebra  with  the  Banach  space  (question  (3)).  The  pullback  algebra,  formed  via  the  topological 
conjugacy,  also  has  spectral  expansions  since  eigenfunctions  are  preserved  when  composing  with  topological 
conjugacies  (question  (2)). 

4)  Proposed  a  novel  operator-theoretic  framework  to  study  global  stability  of  nonlinear  systems. 

Based  on  the  spectral  properties  of  the  so-called  Koopman  operator,  our  approach  can  be  regarded  as  a 
natural  extension  of  classic  linear  stability  analysis  to  nonlinear  systems.  The  main  results  establish  the 
(necessary  and  sufficient)  relationship  between  the  existence  of  specific  eigenfunctions  of  the  Koopman 
operator  and  the  global  stability  property  of  hyperbolic  fixed  points  and  limit  cycles.  These  results  are 
complemented  with  numerical  methods  which  are  used  to  estimate  the  region  of  attraction  of  the  fixed 
point  or  to  prove  in  a  systematic  way  global  stability  of  the  attractor  within  a  given  region  of  the  state 
space. 

5)  Demonstrated  that  the  Koopman  eigenfunctions  and  eigenvalues  define  a  set  of 

intrinsic  coordinates,  which  serve  as  a  natural  framework  for  fusing  measurements  obtained  from 

heterogeneous  collections  of  sensors  in  systems  governed  by  nonlinear  evolution  laws.  These  measurements  can 

be  nonlinear,  but  must,  in  principle,  be  rich  enough  to  allow  the  state  to  be 

reconstructed.  We  approximate  the  associated  Koopman  operator  using  extended  dynamic  mode 

decomposition,  so  the  method  only  requires  time  series  of  data  for  each  set  of  measurements,  and 

a  single  set  of  “joint”  measurements,  which  are  known  to  correspond  to  the  same  underlying  state. 

We  apply  this  procedure  to  the  FitzHugh-Nagumo  PDE,  and  fuse  measurements  taken  at  a  single 
point  with  principal-component  measurements. 

6)  Studied  extreme  phase  sensitivity  in  systems  with  fractal  isochrons.  Sensitivity  to  initial  conditions  is  usually 
associated  with  chaotic  dynamics  and  strange  attractors.  However, 

even  systems  with  (quasi)perlodlc  dynamics  can  exhibit  It.  In  this  context  we  report  on  the  fractal 
properties  of  the  isochrons  of  some  continuous-time  asymptotically  periodic  systems.  We  define  a  global  measure 
of  phase  sensitivity  that  we  call  the  phase  sensitivity  coefficient  and  show  that  it  is  an  invariant  of  the  system 
related  to  the  capacity  dimension  of  the  isochrons.  Similar  results  are  also  obtained  with  discrete-time  systems.  As 
an  illustration  of  the  framework,  we  compute  the  phase  sensitivity  coefficient  for  popular  models  of  bursting 
neurons,  suggesting  that  some  elliptic  bursting  neurons  are  characterized  by  isochrons  of  high  fractal  dimensions 
and  exhibit  a  very  sensitive  (unreliable)  phase  response. 

7)  Developed  a  systematic  approach  to  zoning  and  model  reduction  for  energy  systems  in  buildings  using 
Koopman  Mode  Analysis.  As  the  scope  of  building  design  and  construction  increases  and  building  systems 
become  more  integrated,the  use  of  building  energy  models  has  become  increasingly  widespread  in  evaluating  and 
predicting  building  performance.  Despite  the  growing  sophistication  of  building  modeling  tools,  errors  can  arise  from 
approximations  that  are  made  by  a  practitioner  during  model  creation.  We  examined  the  process  of  model  zoning,  i. 
e.,  how  the  volume  of  a  building  is  divided  into  regions  where  properties  are  assumed  to  be  uniform.  Zoning  is 
performed  during  model  creation  to  decrease  model  complexity.  However,  accuracy  reduces  when  dissimilar 
regions  of  a  building  are  defined  by  a  single  zone.  In  this  paper,  a  systematic  approach  to  creating  zoning 
approximations  is  introduced.  Utilizing  the  Koopman  operator,  the  time-series  output  produced  by  a  building 
simulation  can  be  decomposed  into  spatial  modes  which  capture  the  thermal  behavior  of  a  building  at  different 
time-scales.  Identification  of  spatial  structures  within  these  modes  forms  a  framework  for  the  creation  of  simplified 
models  of  varying  levels  of  granularity.  In  this  paper,  a  detailed  model  is  analyzed,  and  model  accuracy  is  studied 
as  coarser  building  representations  are  created  using  the  introduced  method. 

8)  Applied  the  concepts  abve  to  substantially  reduce  energy  use  in  a  building  on  UCSB  campus.  In  this  work,  an 
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energy  audit  case  study  of  a  LEED  silver  university  building  is  performed  by  decomposing  building  data  into  spatial- 
temporal  modes  known  as  Koopman  modes.  Koopman  modes  are  influenced  by  heat  loads,  from  internal  and 
external  sources,  that  a  building  is  subject  to,  and  analysis  of  these  modes  provides  greater  understanding  of  a 
building's  thermal  behavior.  In  the  presented  application  of  this  framework,  examples  of  energy  waste  are  identified 
including  faulty  equipment  operation,  unnecessary  equipment  usage,  and  HVAC  operating  conditions  requiring  high 
energy  to  maintain.  Through  addressing  the  issues  discovered,  energy  usage  of  the  case  study  building  is  reduced 
by 

13%  with  no  side  effects  seen  in  building  comfort.  Koopman  mode  analysis  has  previously  been  used  to  study  the 
thermal  behavior  of  building  models,  but  has  seen  only  limited  implementation  using  actual  building  data. 

9)  Applied  stability  analysis  using  Koopman  operator  theory  to  power  grid  stability  issues.We  analyzed  nonlinear 
dynamics  and  stability  in  a  network  of  fixed-speed  induction  generators. 

Training  Opportunities:  Students  and  postdoctoral  fellows  participating  in  this  work  have  unique  opportunity  to  do 
interdiscliplinary 

work  between  mathematics  of  dynamical  systems  and  practical  applications  in  energy  systems. 
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isochrons  are  sets  of  points  that  partition  the  basin  of  attraction  of  a  limit  cycle  according  to  the  asymptotic 
behavior  of  the  trajectories.  The  computation  of  global  isochrons  (in  the  whole  basin  of  attraction)  is  however 
difficult,  and  the  existing  methods  are  inefficient  in  high-dimensional  spaces.  In  this  context,  we  present  a  novel 
(forward  integration)  algorithm  for  computing  the  global  isochrons  of  high-dimensional  dynamics,  which  is  based 
on  the  notion  of  Fourier  time  averages  evaluated  along  the  trajectories.  Such  Fourier  averages  in  fact  produce 
eigenfunctions  of  the  Koopman  semigroup  associated  with  the  system,  and  isochrons  are  obtained  as  level  sets  of 
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Abstract:  A  majority  of  methods  from  dynamical  systems  analysis,  especially  those  in  applied  settings,  rely  on 
Poincar  ^e’s  geometric  picture  that  focuses  on  “dynamics  of  states.”  While  this  picture  has  fueled  our  field  for  a 
century,  it  has  shown  difficulties  in  handling  high-dimensional,  ill-described,  and  uncertain  systems,  which  are 
more  and  more  common  in  engineered  systems  design  and  analysis  of  “big  data”  measurements.  This  overview 
article  presents  an  alternative  framework  for  dynamical  systems,  based  on  the  “dynamics  of  observables”  picture. 
We  present  an  overview  of  several  approaches  to  studying  dynamical  systems  using  the  Koopman  operator, 
which  holds  promise  to  resolve  these  issues.  The  dynamics  are  analyzed  by  looking  at  evolutions  of  functions  on 
the  state  space,  rather  than  directly  at  state  space  trajectories.  The  evolution  can  be  understood  by  expanding  the 
function  into  a  basis  of  eigenfunctions  of  the  Koopman  operator.  The  first  approach  is  based  on  the  Koopman 
modes  (KMs),  wh 
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notoriously  difficult,  especially  for  bursting  neuron  models.  In  [W.  E.  Sherwood  and  J.  Guckenheimer,  SIAM  J. 
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